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Sympathetic laser cooling of ions stored within a linear-geometry, radio frequency, electric- 
quadrupole trap has been investigated using computational and theoretical techniques. The simula- 
tion, which allows 5 sample ions to interact with 35 laser-cooled atomic ions, revealed an instability 
heating mechanism, which can prevent ions below a certain critical mass from being sympatheti- 
cally cooled. This critical mass can however be varied by changing the trapping field parameters 
thus allowing ions with a very large range of masses to be sympathetically cooled using a single ion 
species. A theoretical explanation of this instability heating mechanism is presented which predicts 
that the cooling-heating boundary in trapping parameter space is a line of constant q u (ion trap 
stability coefficient), a result supported by the computational results. The threshold value of q u 
depends on the masses of the interacting ions. A functional form of this dependence is given. 
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I. INTRODUCTION 

Laser cooling (!]] has allowed atomic physicists to 
study gas-phase atomic systems at ultra-low temper- 
atures, largely removing thermal effects and allowing 
many atomic systems to be studied at their most fun- 
damental quantum level (e.g. atom interferometry ^|). 
However, although a range of neutral and ionic atomic 
systems has been cooled to sub-Kelvin temperatures with 
direct laser cooling Q, there is a limitation of this ap- 
proach. Conventional laser cooling techniques, both 
Dopplcr H and sub-Doppler require isolated two- or 
three-level systems. If a system has additional low-energy 
meta-stable states or multiple nearly degenerate ground 
states (e.g. rotational structure in molecular ions) then 
the laser-cooling cycle terminates when the system falls 
into a non-laser-cooled level. The solution to this prob- 
lem generally requires one additional re-pumping laser 
system for each meta-stable state (e.g. Ref. ||). Cooling 
of molecules is especially challenging due to their large 
number of accessible rotational states. Understandably, 
there only exist a small number of very complicated theo- 
retical proposals (e.g. Ref. ||) to directly cool molecular 
ions. Therefore, many atomic systems such as those lack- 
ing the necessary energy level structure and virtually all 
molecular systems cannot be directly laser cooled. How- 
ever, a much wider range of atomic and molecular sys- 
tems can be taken to very low temperatures through a 
process known as sympathetic cooling. 

In this technique, a collection of two interacting 
species, the sample species (which can not be directly 
laser cooled) and the laser-cooled species, are confined 
to the same region of space and isolated from their sur- 
roundings through the use of neutral particle or ion traps. 
As the laser-cooled species interacts with the light field 
its temperature falls. When these cold atoms interact 
or collide with the sample species, they extract energy 
from the sample species thus sympathetically reducing 
the sample species temperature. Recently, a variation 



of this approach has been applied to neutrals whereby 
evaporatively-cooled atoms sympathetically cool a sam- 
ple species to produce Quantum Degenerate Gases from 
species whose self-interaction (the interaction between 
one atom with another identical atom) is too weak to 
permit efficient evaporative cooling on its own [Q. This 
type of sympathetic cooling in neutrals depends funda- 
mentally on the short range interactions (collisions) be- 
tween the directly-cooled and sample species, and thus 
the effectiveness of sympathetic cooling depends funda- 
mentally on the particular pair of species selected and 
their associated inter-atomic or atom-molecule interac- 
tion potentials. This makes it rather difficult to produce 
simple and broadly general rules regarding the effective- 
ness of sympathetic laser cooling that would apply to all 
neutral atomic and molecular species. However, when 
considering sympathetic cooling between ions, the domi- 
nant interaction force is Coulomb interactions, which are 
both long range and quite strong. This means that ions 
do not tend to approach each other close enough that 
the short-range interaction potentials, which differ from 
atom to atom and molecule to molecule, become impor- 
tant. Therefore, the effectiveness of sympathetic cooling 
between charged particles tends to be a function of only 
the mass and charge of each of the two ions, allowing for 
rather generalised rules regarding sympathetic cooling of 
charged particles to be deduced. 

This work marks the beginning of a theoretical and 
experimental program to study the processes associated 
with the sympathetic cooling of ions. The longer term 
goals are to develop the tools and expertise to produce 
low-temperature, gas-phase samples in as wide a range 
of atomic and molecular species as is physically possible 
for use in such areas as molecular spectroscopy, spec- 
troscopy of stable and unstable atomic isotopes, and the 
physical implementation of Quantum Information Sci- 
ence protocols. The following pages outline a compu- 
tational/theoretical study of the apparently simple ques- 
tion: What is the range of ion masses that a particular 



1 



laser cooled ionic species can sympathetically cool in a 
quadrupole ion trap? 

Sympathetic cooling of molecular ions by laser-cooled 
atomic ions within a linear Paul trap has been exper- 
imentally demonstrated |^] using ions of similar mass; 
reaching translational temperatures below 100 mK. Fur- 
ther experimental and molecular dynamics (MD) simu- 
lations have shown axial translational temperatures as 
low as 10 mK using a single ion species to sympatheti- 
cally cool over a wider mass range §, but no particular 
or general limits on this range were determined. More 
recently, work related to a novel ion trap in-situ mass 
spectrometry technique under development by the Baba 
and Waki group |l(| has lead to an MD simulation study, 
which illustrated that for particular trap conditions there 
is a lower bound on the sample ion mass (m s ) that can 
be sympathetically cooled. This limit is 0.54 times the 
mass of the laser cooled ion (m c ) JTl| ]. However, no up- 
per bound on the masses that could be sympathetically 
cooled was examined. This work investigates both the 
cooling of higher mass ions, indicating a lack of a con- 
crete upper bound, and examines the question of sym- 
pathetic cooling of light sample ions in a more general 
setting. The latter point lead to techniques to modify 
the cooling threshold value, and provided a general the- 
oretical framework for the heating processes in terms of 
the concept of instability heating. 

II. SYMPATHETIC COOLING OF TRAPPED 
IONS: THE BASICS 

As mentioned in the previous section, sympathetic 
cooling of ions operates through Coulomb interactions 
with laser cooled ions. Theoretically, such interactions 
between ions that are completely isolated from their sur- 
roundings, other than via the laser-cooling beam, will 
eventually bring all of the ions into thermal equilibrium 
at the laser-cooled ion temperature. However, trapped 
ions are not completely isolated because the electric fields 
(and magnetic fields for Penning Traps) used to suspend 
the ions in space can themselves produce a number of 
heating mechanisms (e.g. rf- heating [Q) and it is the 
balance of laser cooling and trap heating that leads to an 
equilibrium temperature. Therefore, depending on the 
relative efficiency of energy extraction through collisions 
with laser cooled ions and trap related heating effects a 
given ion can either increase or decrease its energy. It is 
this balance that leads to the heating-cooling threshold 
observed in Ref. fll|| . In this work, we will concentrate 
on ions stored in a linear-geometry, rf quadrupole-electric 
field trap (commonly referred to as a linear Paul trap) . In 
this type of trap, the ions are contained in the z-direction 
by a DC trapping field, while a rf-oscillating saddle- 
shaped potential traps the ions in the xy-plane p3[ . 

The stability of a single charged particle in an oscillat- 
ing electric quadrupole field is discussed in many refer- 
ences Jl^ |l5|. Highlighting the relevant points, we start 



with a single ion in a linear Paul Trap, where the confin- 
ing potential along the z-axis of the quadrupole is approx- 
imated as harmonic. The ion is subject to the following 
potential in x and y coordinates 

Urf (y 2 - x 2 ) cos(ta) 

®{%,y) = ' W 

where U r f and Q are the amplitude and angular fre- 
quency of the applied rf-electric potential, respectively, 
and ro is the distance from trap centre to one of the elec- 
trodes. Given that F = — eV$, the equations of motion 
for an ion of mass m and charge e in the x-y plane are 
given by 

x + (e/mr^Urf cos(ttt)x = 0, 

y - {e/mrl)U rf cos(fJ% = 0. (2) 

Recasting these equations in unitlcss parameter form 
(£ = Q.t/2, q u = q x = -q y = 2eXJ r $ j m£t 2 r\) we get 

d 2 u 

^-[2^008(20^ = 0, (3) 

where u represents either x or y. Equation (|^) is the 
Mathieu equation in its canonical form Jl6{ , its general 
solution being 

oo oo 

u = dd* C 2n e 2m t+a"e-^ £ C 2n e- 2m ^ (4) 

n— — oo n— — oo 

where a' and a" are integration constants determined by 
the initial conditions, (i.e. u(£o)j w(£o)j an d £o); while 
Ci n and depend solely on the value of q u and not the 
initial conditions Jl4[ . If /i is purely imaginary then the 
ion will have an oscillatory motion (trapped particle). 
However, if /i is real or complex the motion will be un- 
bound (unstable particle) . The determination of /i from 
q u in general is a difficult problem and is best handled 
numerically. It is well known that if \q u \ < 0.908, then /j, 
is purely imaginary and if \q u \ > 0.908 /i becomes com- 
plex, other than for some very localised values of q u [ p"3"| . 
For this reason q u is called the stability parameter. 

For stable trapping parameters, the motion of an ion 
in x and y, as described by Eq. (|4|), can be broken down 
into a superposition of the ion's micromotion and its sec- 
ular motion. The micromotion is an oscillatory motion at 
the angular frequency of the trapping field that increases 
in amplitude the further the ion is from the trap centre. 
It is the direct result of the oscillation of the trapping 
field driving the ion back and forth as the field reverses. 
The lower frequency secular motion is a simple harmonic 
oscillation in the time-averaged pseudo-potential gener- 
ated by the trapping field. The secular frequency of this 
motion, for small values of q Ul is given by [|4| 

^SCC = ""7^^! (5) 
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while the amplitude can be expressed in terms of the 
temperature based on the assumption of ^kgT of energy 
in the relevant mode: 



u 



4 
q u fl 



in 



(6) 



Sympathetic cooling requires that multiple interacting 
ions are simultaneously stored in the trap. We cannot 
deal with such a many-body system analytically, but 
rather it is possible to calculate the motion of each ion 
by integrating the equations of motion for each ion in the 
trap. Additional terms have to be added to Eq. (j|) to 
account for the Coulomb repulsion between the ions, as 
well as the laser cooling [ jL7| . The equation of motion 
for z is similar to the x and y equations, although it is 
simpler because the electric field in this direction is DC. 
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FIG. 1. Temporal plots of the average temperature (see 
text) of 8 amu sample ions and 24 amu laser-cooled ions stored 
in a linear Paul trap calculated under three different trap- 
ping conditions. Line (a) is the sample ion temperatures for 
U rf = 50 V, n = 2vr(2.5 MHz), and q u = 0.71; line (b) is 
the sample ion temperatures for U r f = 100 V, Q = 27r(4.5 
MHz), and q u — 0.53; line (c) is the sample ion temperatures 
for U rf = 30 V, n = 2tt(2.75 MHz), and q u = 0.35. Lines (d), 
(e), and (f) are the laser-cooled ion temperatures correspond- 
ing to the conditions for (a), (b), and (c), respectively. 

Our simulation numerically integrated the motion of 
35 laser cooled Mg + ions interacting with 5 sample ions 
using the Euler-Picard predictor-corrector method with 
0.1 to 1 ns time-steps. The ions are confined radially by 
the potential given by Eq. (|lj) and axially by a harmonic 
well with a force constant of 1 mV/mm. The initial con- 
ditions for the simulations started the ions in random 
positions with the sample ions at a temperature of 40 
K and the laser cooled ions at a temperature of 15 K. 
The laser intensity was selected so that in the absence 



of sample ions, rf-heating balanced the laser cooling and 
the laser-cooled ion temperature would remain the same. 
In addition, the laser intensity was maintained at a low 
enough level to prevent the formation of Wigner or ion 
crystals. These conditions were chosen to be equivalent 
to those used in Ref. [0. The simulation output pro- 
vided the average sample ion and laser-cooled ion tem- 
peratures as a function of time. The temperature was 
calculated by averaging the energy of each ion over 4 cy- 
cles of secular motion and then calculating the average 
ion energy for the 5 sample ions and the 35 laser cooled 
ions, independently. The ion temperatures were deter- 
mined from these average energies using pd| 



(E) = -k B T. 



(7) 



Integration was generally carried out for 30 ms, which 
was usually sufficient to determine whether the temper- 
ature of the sample ions was increasing or decreasing 
with time. First, the laser-cooled ions were interacted 
with a range of more massive sample ions (up to 8 times 
the mass of the laser-cooled ions). This confirmed that 
sample ions heavier than the laser-cooled ions generally 
cooled. However, as the mass difference increased, the 
rate of cooling decreased. 




FIG. 2. Temporal plots of the average temperature (see 
text) of 12 amu sample ions and 24 amu laser-cooled ions 
stored in a linear Paul trap calculated under three different 
trapping conditions. Line (a) is the sample ion temperatures 
for Urf = 40 V, n = 2tt(2 MHz), and q u = 0.59; line (b) is the 
sample ion temperatures for U r f = 70 V, Q = 27r(2.6 MHz), 
and q u = 0.61; line (c) is the sample ion temperatures for 
Urf = 40 V, n = 2vr(2.5 MHz), and q u = 0.35. Lines (d), (e), 
and (f) are the laser-cooled ion temperatures corresponding 
to the conditions for (a), (b), and (c), respectively. 
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Next the case in which the sample ions are less mas- 
sive than the laser-cooled ions was examined. For a given 
U r f and f2, heating or cooling of the sample ion was ob- 
served to depend on its mass. We confirmed that for the 
same trap conditions as in Ref. [O], there is the same 
lower bound on the sample ion mass that can be sympa- 
thetically cooled, i.e. m s — 0.54m c . At this point, tests 
were made to determine if sample ions with m s < 0.54m c 
could be cooled by changing the trap parameters. It was 
found that there always exist values of q u for which this 
is possible. A large number of simulations for a range of 
sample ions with mass below 0.54m c were carried out. In 
all cases the results were more or less the same as those 
illustrated in Figs. [j] andj| Figure |l| is for (m s — 8 amu, 
m c = 24 amu) and Fig. |2]is for (m s = 12 amu, m c = 24 
amu) . Both of these figures show the typical result that if 
the value of q u is reduced sufficiently, the sample ions can 
be cooled. To further test this hypothesis, maps of the 
U r f-VL space were made to determine the trapping field 
parameters for which the sample ions heated, cooled, or 
remained basically unchanged in temperature. Figure || 
shows a pair of typical maps. The q u — 0.908 curve indi- 
cates the single sample ion trapping stability threshold. 
As can be seen, sympathetic cooling does not occur close 
to this threshold and thus the ions are heated. As we 
move away from this threshold we enter a region where 
there appears to be a balance thus the ions neither heat 
nor cool. Finally, for low enough values of q u sympa- 
thetic cooling begins to occur. This begs the question of 
why does cooling occur away from the single ion stability 
threshold and can we predict the threshold for sympa- 
thetic cooling? The answer to these questions are given 
in the next two sections. 



III. INSTABILITY HEATING THEORY 

Although we have not yet discussed the precise heating 
mechanism responsible, the results summarised in Fig. y 
clearly suggest the existence of a second mass thresh- 
old for the behaviour of trapped ions. The first mass 
threshold, dictated by the \q u \ < 0.908 trap stability cri- 
terion, determines if the ion stays in the trap. The sec- 
ond higher-mass threshold, which is also a function of the 
trapping field parameters, determines whether or not the 
sample ion can be sympathetically cooled. 

From the results of the previous section, it is reason- 
able to conclude that the two mass thresholds have simi- 
lar functional dependences on the trapping field parame- 
ters and ion mass, suggesting that the physical causes of 
these two thresholds might be related. These considera- 
tions prompted us to investigate if the heating process is 
related to transient instabilities in the ion motion due to 
interactions with the other ions in the trap, an effect not 
considered in the single ion treatment. We will refer to 
this process as instability heating. 
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FIG. 3. Maps of the regions of heating and cooling in 
U r f — space for 5 sample ions and 35 laser-cooled ions in a 
linear-geometry Paul trap. A indicates simulations that pro- 
duced ion heating, y indicates ion cooling, and • indicates no 
definitive change in temperature. Plot (a) is for 8 amu sam- 
ple ions and 24 amu laser cooled ions while plot (b) is for 12 
amu sample ions and 24 amu laser cooled ions. The solid line 
indicates the q u = 0.908 stability threshold for single ion trap- 
ping, while the second constant q u line in each figure is the 
theoretical threshold between regions of sympathetic cooling 
and heating of the sample ions (see text). 

The basic physical mechanism for instability heating is 
as follows. Equation (|J) gives the general form for mo- 
tion of an ion in an oscillating quadrupole field. For a 
single trapped ion, q u is a constant and for \q u \ < 0.908, 
[i is purely imaginary thus the orbit is oscillatory. When 
\q u \ > 0.908, /i will contain a real part which will cause 
the trajectory of the ion to grow exponentially. However, 
when the ion in question is in the presence of other ions, 
q u can no longer be considered a constant but the addi- 
tional randomly varying Coulomb forces require it to be 
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replaced with a time- varying effective stability parameter 
?eff,«- Viewed as such, |g e ff,u| can now exceed the thresh- 
old value of 0.908 for short periods of times without the 
ion escaping the trap. These time intervals in which the 
ion trajectory grows exponentially are followed by inter- 
vals for which |g c ff,u| drops below 0.908 thus re-stabilizing 
the ion at a larger maximum radius and a higher energy 
orbit. 

To calculate q c s,u, we will consider our physical system 
at a particular point in time and replace the sample ion 
of interest, which is exposed to the trap forces (-Fr>) 
and Coulomb force due to other ions present in 

the trap, with an equivalent ion. Furthermore, we will 
assume that this equivalent ion is positioned at the same 
point in the trap as the ion of interest, having the same 
velocity and acceleration, and is subject to the same trap 
force but is not exposed to the Coulomb force due to 
the other ions, thus requiring it to have a time-varying 
effective mass. The effective mass of the equivalent ion 
in terms of the mass of the sample ion is given by 



rries- 



1 



Ft.u J 



(8) 



Since the effective mass is the mass for an equivalent 
single ion in a trap, Eq. (|^) can be used to obtain g e ff,u 
given by 



where 



q c «,u = <7« (1 + e«) , 



F 

± l,U 



F T ,u 



(9) 



(10) 



Note that g e ff,« is both time varying and can have differ- 
ent values for the x and y directions. We can now write 
the instantaneous stability criterion as 



koff.nl = |g„(l + e„)l < 0.908. 



(11) 



When |<7eff,u| > 0.908 (i.e. when the inter-ionic forces 
sufficiently exceed the trapping forces) we would expect 
the possibility of exponential growth of ion trajectory, 
leading to increased ion temperatures. Figures |]a and 
^b illustrate that this is indeed the case. These simula- 
tions show a short snapshot of the temporal evolution of 
the x-coordinate of a sample ion in parallel with the com- 
puted value of q e s,x- They were obtained with a single 
sample ion, a single laser-cooled ion, and with q u (= 0.9) 
set very near the single ion instability threshold in order 
to maximize the frequency of instability points. The evo- 
lution of the particle motion shows a sequence of steps in 
the amplitude of the ion's motion, corresponding to steps 
in the temperature of the ion. For each of these steps, 
we observe a corresponding spike in q c s yX , where |(? c ff,x| 
exceeds 0.908. In some cases, the amplitude steps down, 
which is expected since Eq. (^) contains both exponen- 
tially growing and decaying terms. The energy drops only 



when the decaying term is much larger than the growth 
term, while it increases under all other conditions. On 
the other hand, there are definitely many points where 
qcS,x exceeds 0.908, but no significant change in the mo- 
tional amplitude occurs. The explanation for this curious 
behaviour is given below. 
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FIG. 4. Temporal plots of the computed evolution of (a) 
q c fi,x (see text), (b) the x-position and (c) the time period that 
q c fi,x remained above 0.908 for the sample ion in a simulation 
involving a 6.6 amu sample ion and a 30 amu laser-cooled ion. 
The trapping field parameters were U r f = 70 V, Q, = 27r(2.9 
MHz), and q u = 0.9. 

What is perhaps more important than q c s jU merely 
exceeding 0.908 is the time interval over which g G ff,u has 
an unstable value. As one can see from Eq. ([Io|), even in 
the presence of a distant laser cooled ion e u undergoes 
a singularity each time the trapping force on the sample 
ion becomes zero. This occurs at least twice per cycle 
of the rf-trapping potential. However, when the time 
period over which q c f{ jU takes an unstable value is short 
compared to the time interval for the rest of the cycle, 
the ion trajectory does not significantly grow, hence no 
measureable heating occurs. Since the ion motion evolves 
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at the secular frequency, a singularity lasting for such 
a short time interval should not significantly effect the 
energy of the sample ion or its temperature. 
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FIG. 5. A temporal plot of the computed trap force on a 
moving ion. Point A illustrates the short instability period 
when the inter-ionic force exceeds the trapping force, while 
period B illustrates the longer instability period occurring 
when the secular motion of the ion passes through x — 0. 

More significantly, e x (or e y ) also undergoes a singu- 
larity as the ion crosses the y (or x) axis. This occurs 
for each half cycle of secular motion. The important 
difference between this and the trapping field singular- 
ities is that the instability periods can last over a much 
longer time period as illustrated in Fig. [|. Here, the solid 
curve represents the trap force on a moving ion as a func- 
tion of time and the horizontal dotted lines represent the 
magnitude of the slowly varying Coulomb force due to 
another ion. For times less than 0.8 /xs, the amplitude 
of the trap force is large (i.e. the sample ion is not near 
the y axis) and thus the inter-ionic force dominates only 
for very short time-intervals as illustrated by A, which 
is too short to significantly effect the trajectory. How- 
ever, when the ion is close to the y axis (i.e. times near 1 
/is) , the trapping force is weaker and the inter- ionic force 
dominates over a longer time period, possibly over sev- 
eral cycles of the trapping field if the ion spacing is small. 
This is illustrated by event B. This results in a measure- 
able heating of the sample ion. The importance of the 
period over which the ion remains unstable is shown in 
part (c) of Fig. ^. We see that the quantized steps in the 
amplitude of the ion motion occur when q e g tX exceeds 
0.908 and does so for a period exceeding 30 to 40 ns. 
To illustrate this point, three events are highlighted. At 
point A, q e fi yX spikes to well above 1 for well over 100 ns, 
resulting in a clear step in the motional amplitude. At 
point B, there is also an amplitude step, but here q e s tX 
barely exceeds 1. However, q e s x stays above the stabil- 



ity threshold for a relatively long time (> 100 ns). On 
the other hand, at point C there are several spikes all ex- 
ceeding 1. However, none of the spikes last longer than 
20 ns, and the ion motion shows no significant change. 



IV. CALCULATION OF HEATING-COOLING 
THRESHOLD 

To generate a predicted shape for the cooling-heating 
threshold, we need to calculate a balance between the 
heating and cooling rates. To calculate the heating rate, 
we first need the variation of the mean instability time 
with the trapping field parameters. From Eq. ( |l(i| ) and 
|TT| , the time interval for which |g e ff,«| is larger than 0.908 
is given by the time interval when 



IF- 



T,u 



< 



(*!,» 



0.908 ± q u 



(12) 



is satisfied. The upper sign is used when Fr tU and F^ u 
have opposite sign and the lower sign when they are 
pointed in the same direction. (F^ u ) is the average inter- 
ionic force during a collision. An estimate of this time 
interval can be obtained as a function of trapping field 
parameters. The trapping force can be determined by 
taking the gradient of Eq. ([!]). Using Eqs. || and [| the 
force on an ion near the trap centre can be approximated 
by (see Fig. |) 

fT,i = — U r f Sin Sit = 2 UrfSmllt , (13) 

where t and t' are used to indicate that the secular motion 
and micromotion may be out of phase. We neglect the 
small effect of micromotion on the x-coordinate of the ion 
near the trap axis. Given that the time scale of interest 
is the secular period, much longer than the trapping field 
period, the trapping field oscillation can be averaged out. 
Strictly speaking this is zero, but stability is dependent 
on |c?eS,ii| < 0.908. Thus it is really the absolute value of 
the forces, not their signs which is of interest. Therefore, 
it is reasonable to replace sin Qt 1 with its rms value. As- 
suming that we are only interested in times close to the 
point where the ion crosses the trap axis, the small angle 
approximation u> soc t -C 1 can be used. Thus 



eu oj soc t 

T> x ~ 9 „2 U rf- 
Lr 



Using Eqs. |5| and || , we find 



Ft,x — 



eU r f /2kBT 
z 'o 



t. 



(14) 



(15) 



Substituting Eq. ( |l5| ) into Eq. (|l2|), and solving for the 
instability time r, we find: 



0.908 2 - ql eU rf V 2fc B T 



0.908|F iiU |g u 4 riK 



(16) 
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One can see immediately that the instability increases 
in duration as q u approaches 0.908, and eventually be- 
comes infinity at the point that q u reaches the single ion 
instability threshold. 

Given Eq. @, we can calculate the heating rate by 
also considering the collision rate, the frequency at which 
these axis-crossing instabilities happen, and the rate at 
which the energy grows during the instability. The col- 
lision rate is simply a function of the number of trapped 
ions and trapping volume, while axis crossings occur 
at twice the secular frequency. To determine the rate 
of energy growth during the instability, we must calcu- 
late how rapidly the radius of the ion trajectory grows. 
The parameter p in Eq. (||) governs this growth. Since 
q u = 0.908 marks the transition from complex to real 
values for p, this coefficient, at least near the transition 
point, should take on the form 



riV\<lu\- 0.908, 



(17) 



where r\ is a constant of order unity. This cannot be 
shown analytically for the Paul trap. However, for the 
closely related rotating saddle trap that can be solved 
analytically, p takes the form of Eq. (|l7|) (l8) . Converting 
back from unitless parameters, the rate of exponential 
growth in the motion during the instability is given by 



1 



1 



p,n=- V nJ\q eS , u \ -0.908, 



(18) 



while the rate of energy increase is simply twice this 
amount due to the harmonic nature of the potential. 
Here, q u has been replaced with q c s, u since it is the quan- 
tity that reaches unstable values. We will approximate 
Eq. (18) using a mean value for g e ff,u by assuming that 
Fi tU is roughly constant over the instability period and 
that Ft, u can be replaced by the average of its maximum 
magnitude during instability and its minimum value of 
zero (sec Eq. (|l2|)). When q c g^ u for parallel and anti- 
parallel forces are averaged, the result is independent of 
q u leaving an average value for the exponential energy 
growth rate of 



(p,Q) = ?7JV0.908. 



(19) 



Finally taking the heating rate to be proportional to the 
the instability time, energy growth rate, the collision rate 
(fcoii), and the axis crossing frequency (2a sec ) yields 



Rh 



K h (X con 0.908\F i!U \q u X] 



0.908 2 



0.908 



Air V m s kBT 



(20) 



To avoid overall heating of the ions, Rh must be bal- 
anced by the cooling effect, which occurs via elastic colli- 
sions between slow and fast moving particles of different 
masses. Simple classical mechanics arguments show that 
the efficiency of elastic collision energy transfer decreases 
as the mass difference between the particles increases. 
The energy transfer efficiency is given by 4p/(l + p) 2 , 



where p = m s /m c |Hj. Taking the cooling rate to be 
proportional to the product of this efficiency and the col- 
lision rate gives 



Re 



m s /m c 



(21) 



When the two rates balance out, we should see a transi- 
tion from heating to cooling of the sample ions. Equating 
Eqs. and EH we obtain 



q„ 



AttK 



0.908 2 



<ll 



0.908|F litl |?7 V 0.908 (1 + m s /m c f 



(22) 



where the constant K, which is independent of the trap- 
ping parameters, is defined to be K c /Kh- The collision 
rates in Eqs. 20 and 21 have been taken to be the same. 
This is valid if the number of laser-cooled ions signif- 
icantly exceeds the number of sample ions. In general, 
cooling occurs when sample ions collide with laser-cooled 
ions, while instability heating is caused by interactions 
between any two ions. However, if most of the ions in 
the trap are laser-cooled ions, then most of the instabil- 
ity heating is a result of interactions between sample ions 
and laser-cooled ions. This means that both the heating 
and cooling rates are proportional to the same collision 
rate, and thus this rate cancels out. Therefore, the rate 
of collisions determines the rate of heating or cooling and 
not the sign of the temperature change. 

Equation ( p2| ) provides the theoretical justification 
that lines of constant q u define the heating-cooling 
boundaries in U r /-il space. Although we cannot calcu- 
late threshold values for q u from first principles, because 
of the unknown scaling factors, it can be done numeri- 
cally. This was done to produce the q u = 0.509 curve in 
Fig. ||b, showing that this line effectively separates the 
regions of heating and cooling. Once the threshold value 
of q u is determined for one mass combination, Eq. ( [22] ) 
can be used to determine the value for other mass com- 
binations. This is illustrated by the q u — 0.408 line in 
Fig. ||a. This line was not fit to the data, but rather it 
was calculated from the constant value for Fig. ||b. It 
should be noted that this model has been applied only to 
a limited mass range and that perhaps other mass effects 
have to be included if very large mass differences are to 
be considered. 

Up to this point, we have concentrated mainly on the 
effect of the mass of the sample ion. However, the laser- 
cooled ion mass is also significant. This is because the 
efficiency of energy transfer between two colliding par- 
ticles decreases as the ratio of the masses deviates from 
unity, becoming effectively zero when one mass is negli- 
gible as compared to the other. However, the instability- 
heating rate is expected to be relatively unchanged as it 
results from the Coulomb interactions between the ions 
and with the trapping field. Thus for larger mass ratios, 
the heating-cooling threshold should move further from 
the single-ion instability threshold in order to reduce the 
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amount of instability heating that must be overcome. 
This reduces the overall collision rate which decreases 
the rate of heating or cooling of the sample ion. The 
effect of laser-cooled ion mass was examined by running 
simulations with 18 amu sample ions interacting with ei- 
ther 24 amu or 36 amu laser cooled ions for a range of q u . 
The results are summarized in Table I which illustrates 
the expected slight shift in the cooling-heating threshold 
for different laser-cooled ion masses. These results are in 
good agreement with the heating-cooling transition val- 
ues predicted by Eq. (23) based on the 12 amu/24 amu 
threshold. The predicted heating-cooling threshold val- 
ues are 0.585 for 18 amu/24 amu and 0.56 for 18 amu/36 
amu. 



V. DISCUSSION AND FUTURE DIRECTIONS 

The results of the preceding sections provide us with a 
physical insight into the competing processes that lead to 
sympathetic heating or cooling of ion mixtures. Gener- 
ally speaking, we can conclude that if sympathetic heat- 
ing is occurring for light sample ions, the trap parameters 
can be adjusted to move away from the single ion insta- 
bility threshold, thereby allowing the sample ions to cool. 
However, there is a limit to this effect. As one moves far 
away from the instability threshold, as would be needed 
to cool very light ions, the trapping field becomes weaker, 
the trapping volume grows, and the collision rate be- 
comes very small. Although an ion still cools, the rate of 
cooling may be so small as to be of little practical use. 

Instability heating is to be distinguished from rf- 
heating p3| , pj[ . This transient effect is caused by 
random interactions between ions, temporarily driving 
them into unstable motion. In contrast, rf- heating of 
thcions near the quadrupole axis is caused by continuous 
Coulomb interactions with the ions which are constantly 
driven by rf field at the periphery. Rf-heating cannot 
be suppressed solely by changing the trap parameters, 
only by reducing the numbers of ions and cooling the 
sample to the point where all the ions crystallize along 
the central axis of the trap where micromotion is not 
present. However, instability heating, which is a major 
source of heating in the low-density gas phase regime, 
can be quenched by choosing favourable trap parameters 
so that the ions are less susceptible to heating. This is 
achieved by minimising the duration of potentially un- 
stable motion and the frequency at which these potential 
instabilities occur, so that in the event of a collision the 
ions is more likely to decrease its energy. 

Further theoretical work is needed to study the effect of 
changing the relative numbers of sample and laser cooled 
ions. When the number of sample ions increases, the in- 
teractions between them can no longer be ignored. This 
will ultimately cause more heating, since such interac- 
tions cannot result in a decrease in the mean energy. Al- 
though, Eq. (16) suggests that instability heating should 



decrease in the higher energy regime this has yet to be 
tested. Such collisions at high temperatures (~300K) 
require much more computational time and other tech- 
niques may be developed to model them. At low energies 
(~1K) Eq. (16), which is applicable for ions in the gas 
phase, becomes less valid since the ions are in the crys- 
talline phase. Instability heating is expected to be in- 
significant here since the axial separation is always large 
enough to prevent the transverse forces, which cause in- 
stabilities, from becoming too large. The dominant heat- 
ing mechanism in this low energy crystalline phase would 
be rf-heating when the number of ions is large. 

In this work, we have ignored the internal degrees of 
freedom of the sample ions, which is reasonable if the 
ions are all atomic, but less so in the case of molecular 
sample ions. The distinction between atomic sample ions 
and molecular sample ions was deliberately overlooked, 
since the focus of this work is the translational degrees 
of freedom. In ignoring the molecular ions internal struc- 
ture, we have a first order approximation as to whether 
or not the sample molecular ions will translationally heat 
or cool. Translational cooling in this approximation is 
seen as a significant component of a complete model of 
the sympathetic cooling. This is because a major energy 
transfer process is expected to be intra-molecular trans- 
fer of energy from rotational and vibrational modes to 
translational degrees of freedom followed by inter-ionic 
translational energy transfer to the laser-cooled atoms, 
as modeled in this work. 



VI. CONCLUSIONS 

The mechanism through which sample ions are heated 
in a linear Paul trap is a direct result of collisions, the 
same collisions that would otherwise bring about thermal 
equilibrium at a lower temperature. A collision can cause 
an ion to become briefly unstable and thus rather than 
energy being transferred to the laser cooled ion for dis- 
sipation in the radiation field, kinetic energy is pumped 
into the ion from the electric fields if its q c s. u is greater 
than 0.908. We have shown that instability heating is 
minimised when q u is small. Since q u is a function of 
mass, trap frequency, and trap depth, although a lower 
mass limit exists below which sample ions will heat, this 
threshold can be changed by modifying the trap param- 
eters. This mass limit corresponds to a second charac- 
teristic threshold value for q u , above which heating will 
occur [jll]. More generally, there exists a three dimen- 
sional space spanned by mass, trap frequency, and trap 
amplitude which is divided into regions of sympathetic 
heating-cooling. Constraining any two of these variables 
will result in an upper or lower limit in the third. 
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TABLE I. Computational observations of sympathetic 
heating or cooling for different trapping field parameters for 
5 sample ions (m c = 18 amu) interacting with 35 laser cooled 
ions (either m c — 24 amu or m c = 36 amu). The sample ion 
behavior was categorized, based on a 60 ms evolution, into 
heating (H), cooling (C), or no observed temperature change 
(N). 
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